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Most  thermoelectric  device  modeling  neglects  the  non-linear  Thomson  effect  in  order  to  develop 
a  closed-form  solution  to  the  governing  heat  equation.  This  simplified  solution  is  beneficial  for  system 
modeling  and  optimization  when  more  intensive  numerical  techniques  are  prohibitive.  An  averaged 
Seebeck  coefficient  is  often  used  in  conjunction  with  the  closed-form  solution  to  incorporate  approxi¬ 
mately  the  effect  of  this  neglected  term  (termed  the  “standard  model”).  While  the  standard  model  has 
been  accepted  in  the  past  for  materials  under  small  temperature  gradients  and  relatively  constant 
Seebeck  coefficient,  there  has  not  been  a  systematic  assessment  of  validity  of  this  modeling  approach, 
especially  for  emerging  materials  and  large  temperature  gradients.  This  work  rigorously  demonstrates 
the  accuracy  and  limitations  of  the  standard  model  through  analytical  derivation  and  comparison  with 
an  efficient  numerical  solution.  It  is  proven  that  the  standard  model  produces  the  exact  module  output 
power  if  an  integral-averaged  Seebeck  coefficient  is  used,  and  also  that  the  standard  model  provides 
a  reasonably-accurate  estimation  of  module  efficiency,  despite  its  limiting  assumptions.  These  findings 
prove  that  the  standard  model  in  fact  incorporates  the  Thomson  effect  with  sufficient  accuracy  that  it 
may  be  used  to  simulate  and  optimize  thermoelectric  systems  as  an  alternative  to  computationally 
expensive  numerical  simulations. 

©  2012  Elsevier  Masson  SAS.  All  rights  reserved. 


1.  Introduction 


Ps  =  IaTE(TH  -  Tc)  -  I2Rte 


(1) 


As  a  nascent  technology,  the  footprint  that  thermoelectrics  will 
have  in  efficient,  sustainable  and  remote  energy  applications  will 
be  determined  by  device-level  cost  and  performance.  To  address 
this  issue,  much  effort  recently  has  gone  into  development  of  high- 
performance  materials  with  little  impact  on  commercially- 
available  modules.  Pivotal  in  the  transition  from  academic 
research  to  commercial  development  is  improvement  of  device¬ 
level  design  tools,  such  as  simple  and  accurate  models,  that  can 
aid  in  the  transition  of  novel  materials  into  practical  power 
production  modules  and  can  also  be  used  for  modeling  of  ther¬ 
moelectrics  in  larger  power  generation  systems. 

The  standard  analytical  model  which  is  the  most  widely  cited  is 
the  seminal  work  published  by  Ioffe  in  1956  [1  ]  and  is  the  basis  for 
other  analytical  models  [2,3].  The  standard  model  (denoted  by 
subscript  ‘s’)  power,  Ps,  and  efficiency,  r]s ,  of  a  thermoelectric  device 
is  typically  in  the  form  [1  ] 
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=  I°lte(Th  -  Tc)  -  I2Rte  , 2) 

S  I^teTh  +  I<te(Th  ~  Tc)  -  \2ReeI2 

where  aEE  is  some  average  Seebeck  coefficient,  I  is  the  device 
current,  KEE  and  REE  are  the  overall  device  thermal  conductance  and 
electrical  resistance,  and  TH  and  Tc  are  the  hot  and  cold  side 
temperatures  across  the  device. 

The  ubiquity  of  the  standard  model  is  most  notable  in  the 
derivation  of  the  dimensionless  figure  of  merit,  ZT,  which  is  the 
most  commonly  used  thermoelectric  performance  metric  and 
defined  as 


ZT  = 


(3) 


where  a  and  l<  are  the  electrical  and  thermal  conductivity  respec¬ 
tively.  Although  the  final  form  of  the  standard  model  is  well-known 
and  it  has  even  been  re-derived  [4],  the  assumptions  made  have 
profound  implications  that  have  been  previously  unexplored.  To 
reiterate,  the  principal  assumptions  are  i)  one-dimensional  trans¬ 
port,  ii)  negligible  contact  resistances,  iii)  material  properties  are 
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constant  over  temperature  and  space  (and,  accordingly,  negligible 
Thomson  effect,  since  the  Thomson  effect  arises  only  from  the 
temperature-dependence  of  the  Seebeck  coefficient).  A  slight 
variation  of  Ioffe’s  standard  model  (Ioffe  assumes  a  constant  See¬ 
beck  coefficient,  olje  =  [«h  +  ad/2,  where  a{T) is  the  temperature- 
dependent  Seebeck  coefficient  and  subscripts  H,  C  and  TE  denote 
hot  side,  cold  side,  and  bulk  terms  respectively)  is  presented  by 
Angrist  [5]  wherein  the  Seebeck  coefficient  is  temperature  aver¬ 
aged  across  the  hot  and  cold  temperatures, 

Th 

«te  =  J  a(T)dT  (4) 

Tc 

where  a  functional  form  of  a  with  respect  to  temperature  is 
assumed  to  be  known.  The  latter  model,  hitherto  referred  to  as  the 
standard  integral-averaged  Seebeck  (SIAS)  model,  is  demonstrably 
more  accurate  than  the  former  in  reporting  module  performance 
for  one  specific  physical  configuration  [6].  As  proven  later  in  this 
paper,  the  higher  accuracy  of  the  SIAS  model  stems  from  the  acci¬ 
dental  incorporation  of  the  Thomson  effect  even  though  the 
authors  neglect  a  term  in  the  governing  equation  that  would 
seemingly  preclude  it. 

Although  the  Thomson  effect  can  be  incorporated  into  a  one¬ 
dimensional  model  via  numerical  techniques,  these  approaches 
are  cumbersome  and  often  not  feasible  in  optimizing  both  indi¬ 
vidual  and  multiple  thermoelectric  module  systems.  Crane  and 
Jackson  assumed  a  constant  Seebeck  coefficient  in  the  optimization 
of  a  cross  flow  heat  recovery  system  [7].  More  recent  high  power 
heat  recovery  system  model  and  optimization  studies  have 
assumed  a  constant  Seebeck  coefficient  [8,9]. 

Some  work  has  attempted  to  account  for  the  Thomson  effect  in 
analysis,  but  only  for  a  limited  case.  Chen  et  al.  proposed  a  model 
that  includes  the  Thomson  effect  and  determine  a  threshold  criteria 
for  neglecting  the  Thomson  effect  based  on  material  properties 
[10].  Assumptions  in  this  model  are  the  same  as  the  SIAS  model 
with  the  exception  that  the  Seebeck  coefficient  has  a  prescribed 
temperature  dependency  (the  other  material  properties  are 
constant).  In  order  to  solve  the  governing  equations,  the  Thomson 
coefficient  is  assumed  to  be  constant,  which  is  a  very  specific  and 
restrictive  case.  The  generalized  governing  equation  for  heat 
transfer  in  a  thermoelectric  material  is 

i  dT  2  /  ri  ft -\ 

kdx2+JXdX+J  /(T  =  0  (5) 

where  k  is  the  thermal  conductivity  and  assumed  to  be  constant  J  is 
the  current  density,  t  is  the  Thomson  coefficient,  and  a  is  the 
electrical  conductivity  and  assumed  to  be  constant.  The  Thomson 
coefficient,  when  held  constant,  leads  to  a  straightforward  closed- 
form  solution  of  the  preceding  equation,  and  furthermore  yields 
a  logarithmic  function  of  temperature  for  a{T).  Freunek  et  al. 
expanded  an  even  more  intricate  model  by  including  contact 
resistances  in  addition  to  the  Thomson  effect,  but  the  model  is  still 
limited  to  a  constant  Thomson  coefficient  in  order  to  solve  the 
governing  equation  [11].  As  such,  in  order  to  maintain  accuracy, 
both  models  are  limited  to  cases  where  the  Seebeck  coefficient  has 
a  logarithmic  functional  form  with  respect  to  temperature,  which 
does  not  correspond  to  many  current  practical  devices,  and  might 
not  be  sufficient  for  emerging  thermoelectric  materials. 

Min  et  al.  experimentally  demonstrated  the  importance  of 
including  the  Thomson  effect  in  the  determination  of  ZT  for  devices 
subject  to  large  temperature  gradients  [12].  Min  et  al.  attempted 
to  account  for  the  Thomson  effect  with  a  simple  model  that  assu¬ 
mes  half  of  the  Thomson  heat  flows  to  the  hot  end  and  the 


other  half  flows  to  the  cold  end.  This  approach  provided  qualitative 
estimates,  but  is  not  sufficient  for  quantitative  analysis.  Yamashita 
examined  the  effects  of  linear  temperature  dependence  on  conver¬ 
sion  efficiencies  [13].  The  approach  assumes  a  linear  temperature 
profile  and  attempts  to  account  for  the  Thomson  effect  using  the 
same  approximation  used  by  Min  et  al.  [12]  that  the  effect  can  be 
accounted  for  by  using  a  Thomson  coefficient  at  the  edge  of  a  ther¬ 
moelectric  element. 

Wee  recently  developed  a  complex  analytical  solution  for 
thermoelectric  devices  with  linear  and  polynomial  material  prop¬ 
erties,  which  accounted  for  the  Thomson  effect  [14].  This  solution 
does  require  the  assumption  of  a  linear  temperature  profile  in  the 
thermoelectric  elements.  Wee  compared  the  analytical  model  to 
numerical  solutions  for  a  half-Heusler  alloy  subjected  to  a  range  of 
large  temperature  differences  with  good  agreement. 

What  follows  is  a  mathematical  account  for  the  Thomson  effect 
and  how  the  SIAS  model  implicitly  incorporates  the  effect  for 
power  predictions  through  the  integral  averaged  Seebeck  coeffi¬ 
cient.  An  efficient  and  accurate  numerical  model  that  solves  the 
nonlinear  thermoelectric  heat  equation  is  used  as  a  benchmark  to 
validate  the  SIAS  model,  Eqs.  (1)  and  (2)  with  Eq.  (4).  Whereas 
power  predictions  of  the  SIAS  model  are  found  to  be  accurate 
enough  for  engineering  applications,  efficiency  predictions  are 
more  sensitive  to  assumptions  in  the  SIAS  model.  For  this  reason, 
a  systematic  parametric  study  relevant  to  potential  new  materials  is 
performed  to  test  the  accuracy  bounds  of  the  SIAS  model. 


2.  An  examination  of  the  Thomson  effect 


In  actuality,  the  Thomson  effect  is  present  in  most  thermo¬ 
electric  configurations.  When  a  thermoelectric  device  with 
a  temperature-dependent  Seebeck  coefficient  experiences 
a  temperature  gradient,  power  is  generated  as  charge  carriers 
respond  to  the  changing  voltage  field  along  its  length  while 
maintaining  a  constant  current.  The  neglect  of  this  non-linear  effect 
without  careful  inspection  of  its  impact  on  thermoelectric  perfor¬ 
mance  can  lead  to  error  in  device  performance  predictions. 
Elowever,  the  influence  of  the  Thomson  effect  is  difficult  to  quantify 
because  the  corresponding  coefficient  causes  the  governing  equa¬ 
tion  to  become  non-linear,  as  seen  in  Eq.  (5).  The  Thomson  coeffi¬ 
cient  is  expressed  as 


T  _  Td«(J) 

dT 


(6) 


For  a  case  where  the  Thomson  coefficient  is  constant  (t  =  Tconst ), 
Eq.  (6)  can  be  rearranged  to  yield 


Ol(T)  —  T Const  In  (r)  +  C 


(7) 


where  C  is  a  constant  of  integration  and  is  left  arbitrary  (Chen  et  al. 
[7]).  To  quantify  the  relative  impact  of  the  Thomson  coefficient, 
a  special  case  is  examined  where  the  logarithmic  functional  form  of 
a(  T)  is  appropriate,  and  thus  t  has  a  single,  constant  value.  Using 
typical  material  property  data  a  BiTe  device,  a  logarithmic  fit  was 
applied  and  the  Thomson  coefficient  was  found  to  have  a  value 
Tconst  =  —0.093  [mVK-1],  as  seen  in  Fig.  1. 

To  compare  the  relative  impact  of  the  Thomson  effect,  the  power 
production  of  a  thermoelectric  material  is  examined.  Chen  et  al. 
[10]  derived  the  power  generation  for  the  limiting  case  of 
a  constant  Thomson  coefficient: 


P Thomson  =  ~  ^CTc)  “  Const  (Th  ~  Tc)  ~  TE  (8) 

where  I  is  the  current  and  Rte  is  the  electrical  resistance  of  the 
thermoelectric  leg.  Considering  the  same  material  with  a  hot  side 
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r/K 

Fig.  1.  Logarithmic  curve  fit  for  Seebeck  coefficient  for  a  typical  BiTe  module. 

and  cold  side  temperature  of  450  K  and  350  K  respectively,  the 
Thomson  effect  (represented  by  the  second  term  in  Eq.  (8))  is  91% 
the  size  of  the  Seebeck  effect  (the  first  term  in  Eq.  (8)).  Given  its 
relative  size,  the  Thomson  effect  plays  a  substantial  role  in  the 
performance  of  a  thermoelectric.  The  difference  between  including 
the  Thomson  effect  and  neglecting  the  second  term  in  Eq.  (8)  can  be 
seen  in  Fig.  2,  where  the  power  of  a  thermo-element  leg  pair  is 
predicted  (assuming  the  material  properties  in  Fig.  1, 
Rte  =  0.0142  Q,  Th  =  450  I<  and  Tc  =  350  I<). 

By  comparison,  the  SIAS  model,  Eqs.  (1)  and  (4),  is  also  shown. 
The  excellent  agreement  between  the  SIAS  model  and  the  Thomson 
effect  model  is  not  accidental.  Combining  Eqs.  (1 ),  (4)  and  (7)  yields 
a  result  identical  to  Eq.  (8). 

Thus,  for  the  case  of  a  constant  Thomson  coefficient  the  power 
predicted  by  the  SIAS  model  matches  the  power  predicted  by  the 
analytical  Thomson  model,  with  a  very  slight  deviation  due  to  the 
error  introduced  by  curve-fitting  the  Seebeck  coefficient  data  with 
the  logarithmic  curve.  This  is  an  unusual  finding  considering  the 
Thomson  effect  was  neglected  in  the  derivation  of  the  SIAS  model;  in 
prior  work,  this  approximation  was  made  without  any  comparison 
to  numerical  models  or  detailed  justification  (see  Ioffe  and  Angrist 
and  others  [1,5,15]).  The  reason  why  the  models  agree  for  power 


production  is  because  when  the  Seebeck  coefficient  is  integrated 
across  the  temperature  field  (as  in  the  case  for  the  SIAS  model),  it 
implicitly  accounts  for  the  Thomson  effect  in  the  single  constant 
term,  ajE •  This  result,  shown  above  for  a  specific  Seebeck  coefficient 
temperature  dependence,  is  actually  general  for  any  temperature 
dependence  as  shown  in  Section  3.  Note  however,  that  while  the 
power  expression  is  precise  with  an  integral  averaged  Seebeck 
coefficient,  the  thermoelectric  efficiency  is  not  (see  Section  5). 

3.  Proof  that  the  standard  integral  averaged  model  exactly 
predicts  power  output 


In  the  previous  section,  it  was  demonstrated  that  power 
predictions  from  the  standard  model  and  the  full  model  are  iden¬ 
tical  for  the  case  of  a  constant  Thomson  coefficient.  It  is  now  proven 
that  this  result  is  valid  regardless  of  the  functional  form  of  the 
Seebeck  coefficient.  Consider  a  single  thermoelectric  leg  of  length  L 
with  constant  thermal  and  electrical  conductivities,  a  current 
density  flowing  along  the  x-axis  of  J,  and  a  hot  side  and  cold  side 
heat  flux  of  qn  and  q'c  respectively.  See  Fig.  3  for  variable  definitions 
of  a  single  leg  pair. 

Eqs.  (5)  and  (6)  can  be  combined  and  integrated  across  the 
length  of  the  thermoelectric  leg  to  yield: 


|»gdx  +  /jr|dx  +  /£**°  (9) 

0  0  0 

Evaluating  the  second  term  using  integration  by  parts  and 
noting  that  J  is  constant,  the  result  is: 

0  \  0 


j(mx=L-lT*]x=o  -j  «dT 


(10) 


\  Tc  J 

After  Eq.  (10)  is  substituted  into  Eq.  (9)  the  result  is  rearranged, 
and  the  first  and  third  terms  in  Eq.  (9)  are  integrated,  the  following 
equation  result  is  obtained: 


kdT 

dx 


-  k 


x=L 


x=0 


+  JP'alx=£.  —  iFalx=0 


(ID 


Qh 


Fig.  2.  Comparison  of  analytical  power  models  including  and  excluding  Thomson  Fig.  3.  Thermoelectric  leg  pair  (thermoelement)  with  relevant  boundary  conditions 

effect  for  constant  Thomson  coefficient  scenario  in  addition  to  the  SIAS  model.  and  sign  conventions.  Heat  transfer  occurs  in  the  x-direction  only. 
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Close  inspection  of  Eq.  (11 )  reveals  that  the  left  hand  side  has  the 
same  form  as  the  expression  for  power  density  of  a  thermoelectric 
material  according  to  the  full  model.  That  is  the  heat  flux  on  the  hot 
and  cold  boundaries  of  a  thermoelectric  leg  material,  as  well  as  the 
power  density  generated  can  be  defined  in  dimensional  form  as: 


„  .  dT 


x=L 


mi  //  1  dT 

®hThJ,  Qc  =  kfa 


x=0 


acTcJ 


(12) 


+  mhThJ,  -®cTcJ  (13) 

x=0 

The  right  hand  side  of  Eq.  (11)  represents  the  power  density 
according  to  the  SIAS  model  (recalling  that  J  =  I/A  and  P  =  P/A)  as 
seen  by  combining  Eqs.  (1)  and  (4).  Simply  put,  the  full  model, 
which  is  based  on  the  non-linear  governing  equation,  and  the  SIAS 
model,  which  ignores  the  non-linear  term  and  compensates  with 
an  integral-averaged  Seebeck  coefficient,  report  the  exact  same 
power: 


p"  =  <&-$  =  * 


-  k 


x=L 


P"  =  P£ 


J  J  adT-J-L  =JaTE(TH-Tc )  -J-L 

Tc 


(14) 


It  is  particularly  remarkable  that  the  expression  for  power 
density  in  Eq.  (14)  does  not  contain  any  reference  to  the  tempera¬ 
ture  gradient.  Regardless  of  the  temperature  field,  the  power 
generated  by  a  thermoelectric  module  can  be  calculated  exactly  as 
long  as  the  functional  form  of  the  Seebeck  coefficient  is  known,  in 
addition  to  the  hot  and  cold  side  temperatures.  The  two  models  still 
predict  different  temperature  fields  (due  to  the  incorporation  or 
neglect  of  the  Thomson  effect),  but  because  no  term  requires 
information  on  T(x),  the  power  generation  predicted  by  the  SIAS 
model  is  exact.  Conversely,  the  calculation  of  thermodynamic  effi¬ 
ciency  does  require  information  from  the  temperature  gradient,  so 
for  such  calculations  the  SIAS  model  will  have  some  error  as  dis¬ 
cussed  in  the  Section  5. 


4.  Dimensionless  modeling 

To  facilitate  comparison  and  validation,  variables  in  the  SIAS 
model  and  the  general  governing  equation  are  made  dimension¬ 
less.  When  Eqs.  (5)  and  (6)  are  combined  with  the  following  non- 
dimensional  expressions: 


the  ratio  of  heat  absorbed  or  liberated  due  to  the  difference  in  the 
Seebeck  coefficient  at  the  surfaces  of  a  leg  pair  that  are  at  Tc  to  the 
heat  conducted  across  a  passive  device  subject  to  a  temperature 
difference  of  AT,  S  represents  the  non-dimensional  temperature 
potential  related  to  the  Carnot  efficiency  (Carnot  Efficiency  =  S/ 
(S  + 1 )),  and  Ra  is  the  ratio  of  the  Joule  heating  to  the  heat  conducted 
across  a  passive  device  subject  to  a  temperature  difference  of  AT. 
The  ratio  6  is  defined  in  Eq.  (18),  but  appears  in  the  non-dimensional 
efficiency  that  follows.  For  the  purposes  of  quantifiably  comparing 
the  full  model  with  the  SIAS  model,  the  non-dimensional  Seebeck 
coefficient  is  prescribed  a  second-order  polynomial  form  as  a  func¬ 
tion  of  temperature  with  arbitrary  constants,  A\,  A2,  and  A3  to 
represent  a  wide  range  of  potential  a  profiles: 

a  =  Ai  +  A2T  +  A3T2  (19) 

For  the  prescribed  Seebeck  coefficient  form,  Eq.  (19),  non- 
dimensional  heat  fluxes  are  derived  using  Eqs.  (18)  and  (12): 


a"  =  fk 

qH  kAT/L 


dT 

dx 


+  £{A\  +  A2  +  A3  +  d)  (S  +  1 ) 


(20) 


d"  - 

Qc  - - 


dT 


/<AT/L  dx 


+  e(A1  +  Q) 


(21) 


x=0 


The  dimensionless  power  flux  and  thermodynamic  efficiency 
can  be  expressed  as  the  following: 


p"  _  q"  q"  _  ^ 

P  -qH~qc-d* 


X=1 


dT 

dx 


x=0 


-\-e(A2  +A3  +S(Ai  -\-A2  +A3  +  6)) 

(22) 


V  =  P"/qfH  (23) 

Eqs.  (22)  and  (23)  provide  the  basis  for  evaluating  the  perfor¬ 
mance  of  a  thermoelectric  and  will  be  used  for  quantitative 
comparison  between  the  models.  A  similar  analysis  can  be  per¬ 
formed  on  the  SIAS  model  whose  non-dimensional  governing 
equation  neglects  the  Thomson  coefficient  (the  ‘s’  subscript  denotes 
the  standard,  or  SIAS,  model): 

d2T 

f2f  +  RA=0  (24) 


=  T  -  Tc  _ 
r  =  T’“  = 


«  -  «Min  ^  _  X 

Aa  ,X  ~  L 


with  boundary  conditions  Ts(x  =  0)  =  0,Ts(x  =  1)  =  1.  The 
temperature  gradient  can  then  be  obtained  by  direct  integration: 


the  dimensionless  general  governing  equation  and  corresponding 
boundary  conditions  (referred  to  henceforth  as  the  full  model) 
becomes: 


(fir 

dx2 


+  e(ST  + 


1) 


dadT 

dTdX 


+  Ra  —  0 


(16) 


1-4  (25) 

The  heat  flux  for  the  SIAS  model  varies  slightly  from  the  full 
model  in  that  the  Seebeck  coefficient  is  assumed  to  be  constant,  as 
expressed  in  Eq.  (4).  The  heat  flux  for  the  standard  model  can  be 
expressed  as: 


T(x  =  0)  =  0,  T(x  =  1)  =  1 


(17) 


where  the  constant  dimensionless  parameters  are  defined  as: 


_  JbaTcL 
kAT  ’ 


J2I? 


akAT 


W-Min 

Aa 


(18) 


+  MteTcJ  (26) 

x=0 

Scaling  the  expressions  with  kAT/L,  and  substituting  the 
temperature  gradient  from  Eq.  (25),  the  non-dimensional  heat 
fluxes  are  expressed  as: 


„  .dT 

qns  ~  kdx 


+  uteJhL  q'cs  =  kfa 


and  AT  =  Th  -  Tc,  A  a  =  aMax  -  &Mim  and  L  is  the  length  of  the  leg 
along  the  x-axis.  aMax  and  aMin  correspond  to  the  maximum  and 
minimum  Seebeck  coefficient  in  the  leg  of  interest.  In  Eq.  (18),  e  is 


q'us  —  1  -  +  £{®-te  +  0)  (5  +  1)?  <7c,s  —  1  +  +  £(aTE  +  0) 

(27) 
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Assuming  the  same  second-order  polynomial  Eq.  (19),  Eqs.  (4), 
(18)  and  (27)  can  be  combined  to  determine  module  power: 


p's  =  - Ra  +  cs(a !  +Y+Y  +  d)  (28) 

Eqs.  (23),  (27)  and  (28)  yield  the  thermodynamic  efficiency,  so 
the  models  have  been  sufficiently  defined  for  comparison. 


5.  Implementation  and  comparison  of  standard  (SIAS)  and 
full  models 


To  solve  the  full  model,  the  non-linear  governing  Eq.  (16)  subject 
to  boundary  conditions  Eq.  (17),  a  numerical  solution  is  created 
using  a  shooting  algorithm.  Eq.  (16)  can  be  split  into  two,  coupled 
first-order  differential  equations: 

-<ST1+1)  %T2~Ra  (29) 


dT 2 
dx 


T2 


(30) 


with  the  starting  initial  value  T^(x  =  0)  =  0.  A  guess  for  the  initial 
value  of  T2  is  made  and  adjusted  using  a  basic  optimization  scheme 
to  ensure  the  problem  boundary  condition,  T \{x  =  1)  =  1,  is 
satisfied.  With  a  good  initial  value  guess  using  the  cold-side  flux 
from  the  standard  model,  only  a  few  iterations  are  required. 
Because  initial  value  problem  algorithms  can  be  used,  large 
matrices  associated  with  finite  difference  and  finite  element  tech¬ 
niques  (methods  typically  used  to  solve  boundary  conditions)  are 
avoided  and  data  storage  is  minimized.  This  drastically  reduces  the 
computational  effort  to  arrive  at  a  solution,  and  retains  high 
accuracy  in  both  the  temperature  field  and  its  gradients,  the  latter 
being  necessary  for  accurate  predictions  of  power  and  efficiency  so 
comparisons  between  the  full  and  SIAS  models  can  be  made. 

The  thermoelectric  leg  pair  power  and  efficiency  for  the  full 
model  and  the  standard  model  are  compared  over  a  range 
dimensionless  parameters:  e  (0.1-1 ),  S  (0.5-2),  RA  (0.1-1 ),  6  (0.1-2) 
and  for  five  different  Seebeck  profiles  A-E  as  depicted  in  Fig.  4. 
Although  the  Seebeck  profiles  do  not  correspond  to  specific  current 
materials,  the  intent  is  to  assess  model  accuracy  for  a  wide  range  of 
potential  profiles  that  could  occur  in  emerging  materials.  The 
ranges  of  dimensionless  parameters  were  determined  to  be 
reasonable  ranges  for  current  and  future  advanced  thermoelectric 
material.  Note  that  with  the  scaling  used,  the  dimensionless  values 
of  a  must  lie  between  zero  and  one  as  indicated. 

However,  as  expected  based  on  the  proof  given  in  section  3,  in 
all  cases  the  power  calculated  for  both  the  full  model  (Eq.  (22))  and 
the  SIAS  (Eq.  (28))  model  were  identical.  The  efficiency  reported  by 
the  two  models  differs  significantly  for  certain  cases.  To  visualize 
where  these  deviations  occur,  the  ratio  77/775  (using  Eqs.  (2)  and  (23)) 
is  mapped  as  a  function  of  e  and  0  for  different  values  of  Ra  and  S. 
The  error  maps  for  Seebeck  profile  A  (Fig.  4)  are  shown  in  Fig.  5(a— 
d)  for  a  range  of  RA  and  S  values.  The  parameters  e  and  6  were 
selected  as  the  independent  variables  because  they  tended  to  have 
the  largest  impact  on  deviations  in  the  models.  As  can  be  seen  in 
Fig.  5(a-d),  the  standard  model  using  the  integral  average  Seebeck 
coefficient  generally  produces  comparable  efficiencies  to  the  full 
model  over  the  ranges  of  dimensionless  parameters  considered. 
The  exception  is  for  the  region  where  e  is  close  to  one,  and  6  is  much 
smaller  than  1.  In  some  cases,  the  error  in  the  SIAS  model  efficiency 
can  be  greater  than  25%  as  seen  in  lower  right  corner  of  Fig.  5(a). 

Similar  results  can  also  be  observed  for  cases  B  and  C,  although 
plots  were  not  included  for  the  purposes  of  brevity.  For  the  Seebeck 


Fig.  4.  Different  cases  for  non-dimensionai  Seebeck  coefficient.  A:  A\  =  0,  A2  =  1, 
A3  =  0.  B:  A-i  =0,A2  =  0,  A3  =  1.  C:  A\  =  0,  A2  =  2,  A3  =  -1.  D:  A\  =0,  A2  =  4,  As  =  -4.  E: 
A\  =  3/4,  A2  =  3/2,  A3  -  -9/4. 


profile  cases  D  and  E,  the  deviation  in  efficiencies  between  the  SIAS 
and  full  models  were  even  less  pronounced  as  shown  in  Fig.  6, 
indicating  the  standard  model  can  be  appropriate  for  estimating 
efficiency  when  the  Seebeck  coefficient  at  the  hot  and  cold  side  of 
the  leg  pair  are  similar  and  have  similar  profiles  to  cases  D  and  E. 

Although  the  errors  in  using  the  SIAS  model  for  profiles  similar 
to  case  A  can  be  significant,  the  entire  region  mapped  in  Fig.  5(a-d) 
are  not  all  physically  possible  nor  of  interest  in  most  applications. 
For  some  cases  the  leg  pair  generation  efficiency  is  less  than  zero, 
which  is  of  little  interest  for  practical  applications.  The  zero  effi¬ 
ciency  is  marked  by  the  white  dashed  line  in  each  figure.  All  points 
to  the  right  and  above  this  line  are  cases  with  efficiencies  above 
zero  (in  the  case  of  Fig.  5(c)  the  white  line  is  not  visible  but  is  in  the 
bottom  left  corner).  Note  that  the  entire  space  in  Fig.  5(c)  has 
positive  efficiencies. 

Current  thermoelectric  materials  have  values  of  ZT  on  the  order 
of  one,  while  some  regions  in  the  figures  are  for  much  larger  ZT 
values.  The  max  ZT  anywhere  in  a  material  with  a  case  A  profile  will 
occur  at  the  hot  side  surface.  Using  the  definition  of  the  dimen¬ 
sionless  parameters,  Eq.  (18),  the  maximum  ZT  can  be  expressed  as 


ZTmax 


e2S(S  +  1)(0+1)2 
Ra 


(31) 


A  ZTmax  equal  to  ten  indicated  by  the  black  bold  line  in  Fig.  5, 
which  represents  a  reasonable  upper  bound  for  the  next  generation 
of  thermoelectric  device  materials.  That  being  said,  a  valid  ther¬ 
moelectric  device  must  fall  between  the  white  dashed  line  and  the 
black  bold  line  in  Fig.  5.  In  some  cases,  the  lines  for  maximum 
Figure  of  Merit  and  zero  efficiency  fall  off  the  figure.  As  can  be  seen, 
generally  the  SIAS  model  is  within  10%  of  the  full  model  for  prac¬ 
tical  and  realistic  thermoelectric  material.  It  is  only  when  6  is  much 
smaller  than  one  and  e  is  greater  than  0.5  that  errors  begin  to  be 
significant,  greater  than  10%,  in  some  cases.  The  largest  errors  occur 
when  both  Ra  and  S  are  large,  as  seen  in  Fig.  5(a).  When  S  is  large, 
there  is  a  high  temperature  gradient  in  the  thermoelectric  leg  and 
any  material  that  has  a  Seebeck  coefficient  that  is  strongly- 
dependent  on  temperature  will  experience  dramatic  changes  in 
the  temperature  field  through  the  thermo-element.  These  dramatic 
changes  cannot  be  accurately  accounted  for  by  the  standard  model, 
so  it  is  expected  that  for  high  values  of  S  there  would  be  greater 
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Fig.  5.  Map  of  efficiency  error  between  SIAS  and  full  model  using  Seebeck  coefficient  dependence  corresponding  to  curve  A  of  Fig.  4.  Filled  contours  represent  the  ratio  of  rj/rjs,  the  black 
line  represents  the  ZTMax  =  10  value  for  the  thermoelectric  and  the  white  dashed  line  represents  zero  efficiency,  (a)  RA  =  1, 5  =  2,  (b)  RA  =  1, 5  =  0.5,  (c)  RA  =  0.1,  S  =  2,  (d)  RA  =  0.1, 5  =  0.5. 


error  in  the  SIAS  model.  Additionally,  with  larger  S  comes  a  greater 
Carnot  efficiency,  or  larger  potential  for  converting  energy. 

Alternatively,  with  larger  Ra  comes  decreased  efficiency,  since 
the  parasitic  Joule  heating  becomes  dominant.  At  smaller  values  of 
S  and  larger  values  of  Ra,  inefficiency  dominates  and  much  of  the 
map  falls  below  the  zero-efficiency  floor  depicted  by  the  white 


Fig.  6.  Worst-case-scenario  map  of  efficiency  error  between  SIAS  and  full  model  using 
Seebeck  coefficient  dependence  corresponding  to  curve  D  of  Fig.  4.  Filled  contours 
represent  the  ratio  of  77/77$,  and  the  white  dashed  line  represents  zero  efficiency.  RA=  1, 
S  =  2. 


dashed  line.  As  such,  the  standard  model  and  the  full  model  have 
fair  agreement,  with  errors  <10%,  as  seen  in  Fig.  5(b). 

As  mentioned  earlier,  the  largest  errors  occur  when  5  and  Ra  are 
both  large.  For  cases  where  an  and  ac  values  are  similar,  such  as  the 
case  Seebeck  profile  case  D,  the  error  is  usually  negligible.  For  the 
worst  case  scenario  using  case  D  less  than  8%  error  is  reported,  as 
seen  in  Fig.  6.  For  most  viable  configurations  in  this  case,  the  error 
would  be  limited  to  2-4%.  All  other  permutations  of  S  and  RA  yield 
even  smaller  errors,  indicating  that  for  this  particular  range  and 
Seebeck  profile  case,  the  standard  model  is  reasonably  accurate. 

Case  E  behaves  slightly  differently  than  case  D  in  that  it  main¬ 
tains  errors  of  less  than  5%  (in  the  worst-case  scenario),  but  rather 
than  over-predicting  efficiency  case  E  under-predicts  it.  The  under¬ 
prediction  likely  stems  from  the  fact  that  case  E  is  the  only  scenario 


Table  1 

List  of  thermoelectric  parameters  similar  to  a  commercial  device. 


Parameter 

Value 

Units 

Aa 

200 

|iVI<-1 

ac 

400 

pVK'1 

a 

60,000 

Sm-1 

k 

1.6 

Wm-1I<-1 

L 

0.001 

m 

Tc 

300 

K 

Th 

450 

K 

J 

1.2  106 

Am-2 

S 

0.5 

- 

e 

0.3 

- 

Ra 

0.1 

- 

6 

2.0 

- 
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where  an  is  less  than  ac,  however  the  relatively  small  error  is  still 
attributable  to  the  fact  that  the  hot  and  cold  side  Seebeck  coeffi¬ 
cients  are  relatively  close  in  value. 

To  demonstrate  the  utility  of  the  error  maps,  parameter  values  are 
chosen  to  correspond  to  those  in  the  range  of  a  typical  thermoelectric 
material  that  is  commercially  available.  Table  1  provides  material 
property  and  module  parameters  as  well  as  associated  dimension¬ 
less  parameters  calculated  from  Eq.  (18).  The  current  density, 
J,  corresponds  to  1.2  A  of  current  given  a  1  mm  by  1  mm  leg  cross 
section.  Since  S  =  0.5  and  Ra  =  0.1,  Fig.  5(d)  is  utilized  to  calculate  the 
error  in  the  model,  and  with  e  =  0.3  and  6  =  2 ,  the  error  in  the  effi¬ 
ciency  calculated  by  the  SIAS  model  is  less  than  2%.  For  this  case, 
then,  the  error  associated  with  using  the  SIAS  model  is  established, 
and  would  be  deemed  accurate  enough  for  engineering  applications. 

6.  Conclusions 

The  Thomson  effect  is  demonstrably  present  in  thermoelectric 
materials  and  often  contributes  significantly  to  the  power  genera¬ 
tion  capabilities  of  a  module,  contrary  to  assumptions  in  modeling 
that  have  been  used  for  decades.  This  work  proves  by  numerical 
study  and  mathematical  derivation  that  the  standard  model  that 
uses  an  integral-averaged  Seebeck  coefficient  (but  neglects  the 
Thompson  term  in  the  temperature  field  evaluation)  incidentally 
calculates  the  power  output  of  a  leg  pair  exactly,  since  the  integral 
average  Seebeck  coefficient  is  consistent  with  the  incorporation  of 
the  Thompson  effect.  Although  the  averaging  of  the  Seebeck  coef¬ 
ficient  and  neglect  of  the  Thomson  effect  is  formally  not  self- 
consistent,  the  SIAS  model  also  estimates  the  efficiency  for  a  ther¬ 
moelectric  reasonably  accurately.  Thus,  the  standard  model  using 
an  integral-average  model  is  a  highly  convenient  tool  that  appears 
to  be  sufficiently  accurate  for  engineering  applications  such  as 
system  modeling  and  optimization,  where  intensive  computation 
simulations  of  individual  thermoelectric  subunits  is  not  practical. 
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Glossary 


Symbol,  Description  and  Unit 

A/,  A2,  A3:  Polynomial  coeff.  of  non-dimensional  Seebeck  coeff.,  — 

A:  Leg  cross-sectional  area,  m2 

I:  Electrical  current,  A 

J:  Electrical  current  density,  Am-2 

k:  Thermal  conductivity,  Wm-1K-1 

I<te:  Thermal  conductance,  WK'1 

I;  Leg  length,  m 

P:  Power  produced  by  thermoelectric  (full  model),  W 

P's :  Power  produced  by  thermoelectric  (SIAS  model),  W 

p":  Power  density  (full  model),  Wm“2 

P"s:  Power  density  (SIAS  model),  Wm'2 

qH:  Hot-side  heat  flux  (full  model),  Wm-2 

qH,s:  Hot-side  heat  flux  (SIAS  model),  Wm-2 

q’c  Cold-side  heat  flux  (full  model),  Wm-2 

qctS:  Cold-side  heat  flux  (SIAS  model),  Wm-2 

Rte:  Electrical  resistance,  Q 

Ra:  Ratio  of  Joule  heat  to  heat  conduction,  — 

S:  Temperature  potential,  — 

T:  Temperature,  K 

Th:  Hot-side  temperature,  K 

Tc:  Cold-side  temperature,  K 

x:  Distance  along  thermoelectric  leg,  m 

ZT;  Figure  of  Merit,  — 

ZTMax-  Maximum  Figure  of  Merit,  — 

Greek  Symbol,  Description  and  Unit 

a(T):  Seebeck  coefficient,  mVK-1 

ac  Hot-side  Seebeck  coeff.,  mVK-1 

a^  Cold-side  Seebeck  coeff.,  mVK-1 

UMax-  Maximum  value  of  Seebeck  coeff.,  mVK-1 

aMin •  Minimum  value  of  Seebeck  coeff.,  mVK-1 

ajc  Integral-average  Seebeck  coeff.,  mVK-1 

77;  Thermodynamic  efficiency  (full  model),  — 

qs:  Thermodynamic  efficiency  (SIAS  model),  — 

e:  Ratio  of  Thomson  heat  to  heat  conduction,  — 

a:  Electrical  conductivity,  Sm-1 

t:  Thomson  coefficient,  mVK-1 

^Comd  Constant  Thomson  coeff.,  mVK-1 

6:  Ratio  of  minimum  to  difference  in  Seebeck  coeff.,  — 


